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Abstract 

Coupled TchebyschefT maps have recently been introduced to explain param- 
eters in the standard model of particle physics, using the stochastic quantisation 
of Parisi and Wu. This paper studies dynamical properties of these maps, finding 
analytic expressions for a number of periodic states and determining their linear 
stability. Numerical evidence is given for nonlinear stability of some of these states, 
and also the presence of exponentially slow dynamics for some ranges of the param- 
eter. These results indicate that a theory of particle physics based on coupled map 
lattices must specify strong physical arguments for any choice of initial conditions, 
and explain how stochastic quantisation is obtained in the many stable parameter 
regions. 

1 Introduction 

Coupled map lattices consist of a continuous variable (here $) defined on a discrete 
lattice representing space and time. The dynamics relates the value of $ at each space- 
time point to its value at previous lattice points, typically just the most recent time and 
the nearest neighbours in space. Coupled map lattices include the simplest examples of 
spatiotemporal chaos and are often used to model less tractable systems such as nonlinear 
partial differential equations. 

Diffusively coupled Tchebyscheff maps were introduced by Beck in order to chaotically 
quantise field theories @, later applied to spontaneous symmetry breaking of quantised 
Higgs fields || . The discrete time of the coupled map lattice corresponds to the fictitious 
time in the Parisi- Wu stochastic quantisation |3[], and is taken to infinity in order to 
reproduce quantum mechanics. The function appearing in the map is naturally given by 
the derivative of the quartic double well potential, which (appropriately scaled) is the 
third degree Tchebyscheff polynomial discussed below. In addition, Tchebyscheff maps 
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are conjugated to Bernoulli shifts, and so have very chaotic statistical properties, hence 
the possibility that the stochastic forces arise intrinsically from the chaotic dynamics. 

In a later work |J, Beck considers coupled map lattices based on the second degree 
Tchebyscheff polynomial in addition to the third degree case considered previously, and 
restricts attention to one spatial dimension. He numerically computes averages of interac- 
tion and self energies, and reproduces to about four digits a large number of parameters 
occurring in the standard model of particle physics, with predictions for some poorly 
constrained quantities such as neutrino and Higgs masses. 

The opinion of the present author is that although Ref. H does not contain a complete 
theory relating particle physics to coupled map lattices, the probability of obtaining such 
a large number of accurate predictions with so few adjustable parameters is sufficiently 
small that the models warrant further attention. 

Quite apart from the context in which these models were introduced, they have a 
number of appealing features from a dynamical point of view. In particular, there is 
an adjustable parameter a in which the system takes two exactly solvable limits, a = 
which is decoupled and fully chaotic, and a = 1 in which the lattice decouples into two 
sub-lattices, the latter of which are fully chaotic for the advanced coupling case and fully 
stable for the backward coupling case (see below for definitions). 

Motivated by the desire for exact results and the simple form of the equations, this 
paper analyses the simplest periodic states, those for which the field does not depend on 
the spatial variable, called "synchronised states". Finding these states is equivalent to 
finding periodic orbits of a one dimensional map, but the linear stability is a more in- 
volved calculation than one dimensional maps due to spatially dependent unstable modes. 
Stability analysis is important because the proliferation of stable states undermines the 
argument that stochastic quantisation arises from the chaotic dynamics. 

A general theory of linear stability of periodic solutions of coupled map lattices is given 
in Refs. f| [/]]. Our analysis of synchronised states leads, as in these works, to cyclic 
matrices, which are always explicitly diagonalisable. The general theory also permits 
spatially periodic solutions of length greater than one (leading to block cyclic matrices), 
and lattices of spatial dimension greater than one (leading to blocks within blocks). 

Finding that an orbit in a spatially extended system is linearly stable does not imply 
that a finite measure of initial conditions approaches the orbit, unlike finite dimensional 
dynamical systems. Actually, there are many possible inequivalent measures of initial 
conditions: we choose an initial distribution that is Lebesgue in the field $ and spatially 
independent, but physical arguments may justify an alternative choice. In fact, in view of 
the results, we will argue that the final (ie infinite fictitious time) state depends strongly 
on this choice of initial conditions, and hence that physical arguments are required to 
specify such a choice. 

The linear stability of the synchronised states is given in Sec. |2|. Numerical work in 
Sec. |3] reveals properties of the dynamics far from these periodic states. The conclusion 
summarises the results, discussing in more detail the implications for the Beck theory. 
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T 3 ($) 
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Table 1: Definitions of functions / and g appearing in Eq. ([!]) for the coupled map lattices 
considered in this paper. The Tchebyscheff polynomials are Tjy($) = cos(iVarccos $), 
specifically T 2 ($) = 2$ 2 - 1 and T 3 ($) = 4$ 3 - 3$. 

2 Linear stability 
2.1 Formalism 

As in low dimensional dynamical systems (such as the individual maps), linear stability 
analysis in coupled map lattices consists of determining whether and how fast infinitesimal 
perturbations of a reference trajectory (here a periodic orbit) grow with time. As noted 
above, however, linear stability, in that all sufficiently small (everywhere) perturbations 
decay with time, does not imply that a finite measure of initial conditions is attracted 
to the reference trajectory. Notwithstanding ths caveats of interpretation, the linear 
stability analysis is a good starting point for finding stable orbits and their corresponding 
parameter ranges; nonlinear stability is discussed in Sec. |3|. 

Consider a diffusively coupled map lattice, in which x and t are integer space and 
time variables respectively (a more transparent notation than the conventional i and n 
respectively) : 

^, t+1 = (1 - a)f($ x , t ) + | [g($ x+1 , t ) + g($ x -i, t )] (1) 

Here, / is the map at each site x, g couples neighboring maps with a strength a. Later 
we will specify / and g to be the possibilities enumerated in Tab. [1]. We demand that 
$ G [—1,1]; this property is preserved if / and g are closed on this interval, and a G [0, 1]. 
Any periodic orbits found by solving the appropriate equations are only relevant if they 
lie in this range. 

Synchronised states are solutions in which <& x t is independent of x, and so are given 
by the one dimensional map 

$ t+1 = (l-a)/($ t ) + a(K*t) (2) 

which is a weighted mean of the / and g maps. To a first approximation, stable solutions 
of the full coupled map lattice which are periodic in time occur when the synchronised 
map (0) has a superstable periodic orbit, that is, when one of its critical points is mapped 
back to itself after some time, leading to a Lyapunov exponent of — oo. We will observe 
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that the linear stability properties of the full coupled map lattice are slightly different 
from that of the synchronised map, so that although stability occurs in the region of 
synchronised superstable orbits, the superstable orbit may not be linearly stable in the 
extended system. 

The difference between stability in the synchronised map and the extended system is 
that the full coupled map lattice has an infinite number of degrees of freedom, leading to 
a richer spectrum of possible instabilities. Differentiating Eq. ([!]) leads to the evolution 
of the perturbations 



or in matrix form 



<5$ m = J t 5$ t 

where Jt is a cyclic (assuming periodic boundary conditions) tridiagonal matrix 
( (1 - a)/'($ t ) Ig'^t) %g'{$ t ) \ 



(3) 
(4) 



Jt 



1 - a)/'(* t ) 

§</($*) (l-a)/'(* t ) 



(5) 



§</(<!>*) (1 - o)/'(*t) j 



The full stability matrix for a periodic orbit of (temporal) length r is a product of the 
form 

J = Jt-\Jt-i ■ ■ ■ Jo (6) 



For a cyclic matrix M, 



'•j 



rrii 



.j of size X (the full extent of the system in the spacelike 



direction), there are X eigenvectors Vk, given by 

v k = (l,e ik ,e i2k ,...,e- ik ) (7) 

where k = 2irj/X, j = 0, 1, . . . , X—l. This, together with the fact that the corresponding 
eigenvalues are 

x-i 



3=0 



follows by substitution into Mvy. = XkVk- 

The product of cyclic matrices J given by 
eigenvalue is 



r-l 



acting on v k thus gives Xf-v k where the 

(9) 



ft [(1- a) f'($ t ) + ag'($ t ) cos k] 
t=o 

Linear stability follows if |Afe| < 1 for all values of k. Note that this is a stronger condition 
than stability of the synchronised map (|2|), which is \Xk\ < 1 for k = 0. 

There are two classes of analytic results obtained using Eq. (^), r < 2 in which case 
the periodic orbits can be given exactly, and certain families of periodic orbits which 
are given exactly in the limit r 
respectively. 



oo. These are investigated in sections |2.2| and ^3 



4 



Dynamics 
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w 
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Table 2: Synchronised states of period r = 1. The first column gives the model as defined 
in Tab. [I]. The second column gives the solution of Eqs. (0,0). The third column gives 
values of a for which $ G [0, 1]. The fourth column gives values of a for which |A^| < lVfc 
in Eq. (^), that is, the orbits are linearly stable. 

2.2 Synchronised states of periods one and two 

The equation for the synchronised dynamics (^) leads directly to states of period r by 
imposing the condition 

$ f+ r = $f (10) 

For r = 2 it might be supposed that such an equation may have very complicated, 
possibly non-closed form solutions, since for the cubic maps the relevant polynomial is 
of degree 9 (albeit reducible to degree 6 by noting that the fixed point r = 1 solutions 
can be removed). It turns out that, even in this case, the solution can be written using 
arithmetic operations together with at most two square-root operations. 

The periodic orbits are given in Tables 0, || We start by noting that maps 2A and 
3v4, for which the synchronised maps @ are just the original Tchebyscheff maps, have 
no stable periodic orbits as these are fully expanding when transformed to the variable 
cf) = arccos$. The fixed points (r = 1) of the 2B and 3B maps are similar in that they 
are always unstable, and in that the solution $ does not depend on a. The remaining 
periodic orbits for r < 2 have regions of stability in a. 

The calculations here are quite straightforward, so details will be given for only the 
most complicated case, the r = 2 orbits of the 3B model. The solutions of the equa- 
tions (0, [10|) are obtained using a symbolic manipulation package, and are given in Tab. |3] 
and plotted in Fig. [I], where the three orbits are designated "central", "upper" and 
"lower". Equation @ reads 

A fe = [3(1 - a) (4$^ - 1) + a cos k) [3(1 - a)(4$ 2 - 1) + acosk] (11) 
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Dynamics $ Domain of validity Domain of stability 

2A \0j] ~ None 

2A~ -i±V5-32a+32q i [ Q , ^ w 0. 1938] (2^5 « 0.1220, 1/6 « 0.1667) 

w 0.8062, 1] (5/6 « 0.8333, w 0.8780) 

2B -i-a±V5-i8a+ 9ai [0,1/3^ 0.3333] (^ ss 0. 1835, « 0.2404) 

2£- ~ 1+a± l 5 r 4 14a+ — [0, 5/9 « 0.5556] ~ 0.3675, 5/9 P3 0.5556) 

3 A ±^ [0,1] None 

3 A ^fl [0,1] None 

3 A -^fi [0,1] None 

3B ±\/fii [0,1/2 = 0.5000] (1/4 = 0.2500,2/5 = 0.4000) 



3B ± ^ 3-4a±V5-24a+16a^ ([0,1/4 = 0.2500] [^p± ^ 0.2118, 1/4 = 0.2500) 

35 ± / 3-4a Tv g4 a+16al [0,1/4 = 0.2500] (^^^ 0.2118,1/4 = 0.2500) 



Table 3: Synchronised states of period r = 2. Columns as in Tab. |2] except that the two 
values of $ are represented by the ± signs. The 2A~ model has two disjoint domains of 
validity of the given orbit. 



where 

l-2a±y / b-2Aa+l&a I 



4$1 - 1 



2=2^T U PP er 

central (12) 



l-2a=FV5-24a+16a 2 l nwpr 
2-2a luwci 

The maximum value of A & is determined from 

= = - a sinJfe[3(l - a) (4$^ - 1 + 4$ 2 - 1) + 2a cos k] (13) 

for which possible solutions are k = 0, k = tt, and 

3(l-a)(4$ 2 -l+4$ 2 -1) ( "3^ upper/lower 

cos A; = = < (14) 

2a -3^ central 

V a 

These solutions can be substituted back into Eq. ([[]]) taking care to exclude ranges of a 
in which cosfc is not in [—1,1]. 

For the central orbit, the k = solution gives 

\ k = 9 _ 48a + 64a 2 (15) 

which has magnitude greater than one for a < 1/4. The k = n solution gives 

A fc = 9 - 60a + 100a 2 (16) 
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a 

Figure 1: Period 2 orbits of the 3B model. The three orbits (see Tab. |3|) are shown using 
solid, dashed and dotted lines for the "central" , "upper" and "lower" orbits respectively; 
in each case there are two values of $. 
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which has magnitude greater than one for a < 1/5 and a > 2/5. The solution cos A; = 
—3(1 — 3a) /a is in the range [—1,1] for a > 3/10 and gives 

X k = (17) 

Thus the domain of linear stability for the central orbit is (1/4, 2/5) as given in Tab. |3|. 
For the upper and lower orbits, the k = solution gives 

A fc = -9 + 48a - 32a 2 (18) 

which has magnitude greater than one for a < (3 — v / 5)/4 ~ 0.1910. It is also equal to 
one at the maximum value of a for the orbits, a = 1/4. The k = n solution gives 

\ k = _9 + 42a - 20a 2 (19) 

which has magnitude greater than one for a < (21 — v281)/20 ~ 0.2118. The other 
solution cosk = —3(1 — 2a) /(2a) is never in the range [—1,1] for a < 1/4, and so is 
irrelevant. Thus the domain of linear stability for the upper and lower orbits is ((21 — 
v^8l)/20, 1/4) as given in Tab. [3]. The linear stability for all orbits is presented in Fig. ^. 

In this manner all the other domains of linear stability in Tables § and |3] were gen- 
erated. Apart from the 2A and 3^4 models for which the synchronised dynamics is 
unaffected by the coupling, there are analytically tractable linearly stable domains of 
period 2 for all models. Among the interesting features are the bifurcation observed in 
the 3B model (see Fig. [I]), and two regions of stability in the 2A_ model (see Tab. |3|). 
We now move to the other analytically tractable case, the limit r — > oo. 

2.3 Synchronised states of long period 

This section gives an analytic treatment of families of synchronised stable orbits with long 
period. The idea is similar to that of the superstable orbits of (for example) the logistic 
map, in which the critical point (at which the derivative is zero) maps back to itself after 
a number of iterations, leading to a periodic orbit with a Lyapunov exponent of — oo. 
Here the extended nature of the dynamics complicates and destabilises the dynamics to 
a degree, however linear stability is still possible, as we shall show. 

In the original Tchebyscheff maps the critical points are pre-images of unstable fixed 
points at $ = ±1, and so never return to form a superstable orbit. However, the synchro- 
nised map (|2|) is, for small a, and excluding the 2A and 3A small perturbation 
on the Tchebyscheff map, so that after a "long" time (of order — In a) the trajectory can 
return to the critical point. There are many possible ways in which this can happen; 
we focus here on the simplest in which the trajectory, after at most a single point near 
$ = — 1 always remains in the rightmost branch of the map. An example of such a 
periodic orbit is given in Fig. 

We now locate a periodic orbit of length r, and in particular calculate the value of a 
at which it is close to superstable. Suppose that for some value of a, the initial point is 
close to the critical point, that is, 

$o = 0(a) (20) 
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Figure 2: Linear stability of the orbits shown in Fig. |l|. Depending on the value of k, 
takes values between the solid or dashed lines for the central or other orbits respectively. 
The orbits are linearly stable if all for a given orbit lie between the dotted lines. 
Formulas for these lines are in Eqs. (P^-19|). 
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Figure 3: Stable r = 4 orbit of the 2B~ synchronised map, Eq. (|2|). In this case a = 
0.0171, slightly greater than the value a* /N T = 0.0165 predicted by the large r theory 
(see Tab. [| below). 
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Model bx b 2 d b* a* n* g 

2A_ 2 10 2 2/3 3tt716 « 1.851 l/(2vr) 

2fi 1 6 3/8 tt 2 /3^ 3.290 1/(2tt) tt/6 « 0.5236 

2fi_ 1 4 2 7/24 3tt 2 /7^ 4.230 1/(2tt) 3tt/14 « 0.6732 

35 3/2 27/2 1/6 4vr 2 /3 « 13.159 1/(4ttv%) ^(3^) « 0.6046 



Table 4: Coefficients of the long period orbits, defined throughout Sec. 2.3 



for the second degree maps (that is, 2A , 2B and 2B ) and 

$o = -1/2 + 0(a) (21) 
for the third degree map 3B. We can use Eq. (H) to find the successive iterates. We have 

$! = -1 + b 1 a + 0(a 2 ) (22) 

for the second degree maps and 

$! = 1 _6 ia + C(a 2 ) (23) 

for the third degree map. Here &i is a known constant that depends on the map, and is 
given in Tab. £|. Iterating again, we have 

$ 2 = 1 - b 2 a + 0{a 2 ) (24) 

for all maps. Again, b 2 is given in Tab. |] for each model. 

Further iterates remain near the unstable fixed point of the unperturbed map at 
$ = 1, that is, we have 

= 1 - bja + 0(a 2 ) (25) 
where each successive bj is obtained from the previous one by an equation of the form 

b j+1 = N 2 bj + d (26) 

where N is the degree of the map (two or three) and d is given in Tab. ^ for each map. The 
recursion (PE| ) is a linear difference equation, solved by multiplying through by N~ 2 ^ +1 > 
and summing over j from 2 to m — 1. The result is 

N 2{m-2) _ i 

b n = jV 2(m ~ 2) 6 2 + — « iV 2m 6* (27) 

where 

6 = A -: ' V ( > _ yi 

is also given in Tab. |j. 



! ' 3 + ^ (» 



11 



As the trajectory approaches the end of the periodic orbit, $ deviates substantially 
from 1, invalidating the above approximation. This part of the trajectory is obtained 
by iterating $o = backwards. Recalling that the Tchebyscheff polynomials can be 
represented as T n ($) = cos n arccos $ we have for the degree 2 models 



7T 



$ T _i = cos - + O(a) 



71 



$ T _2 = cos - + 0(a 



7T 



2 n+l 

and for the degree 3 model 



cos— T + 0(4""a 




+ 0{4T n a) 



2n 



$ T _„ = cos— T + 0(9- n a) 



(29) 
(30) 

(31) 
(32) 



Note that the 0(a) corrections are suppressed by the contraction of the inverse map. 

The orbit is closed by matching the expressions for $ m and $ r _„, assuming that both 
m and n (actually N 2m and N 2n ) are much greater than one. The result is 

1 _ b * 4 r-n a + 0{A -n a) = 1 _1 ' + ^^n) (33) 



for the second degree models and 



1 - b*9 7 



0(9~ n a) 



1 



2tt 
3^+1 



+ 0{9 



-2n\ 



for the third degree model. The value of a can now be read off as 

a = a*N- 2r + 0(N~ At ) 



(34) 



(35) 



where a* is given in Tab. |J The result is consistent with the original assumption that 
a is small. Note that the original 0(a) freedom in the initial condition $ contributes 
to 0(a 2 ) when the trajectory is near $ = 1, and so is expected to shift a by an amount 
of order a 2 , that is, N~ 4t . Obtaining the appropriate coefficient would require a more 
involved calculation. 

Now we turn to the linear stability of these long orbits using Eq. (0). is now a 
product of a large number (r) of factors. The t — factor is small since both /'($o) 
and a are small, and will be considered separately. For the remaining factors the /'(<&) 
term dominates and the g'(<$>) term can be ignored, owing to the presence of the small 
quantity a. Thus we have 



T-l 



A fc = [/'($ ) + ag'(® ) cos k] ]J f($ t )[l + O(ar)} 



(36) 



t=i 
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The product is 



n /'<*.) = - n 4 cos ^ = - n 2 = -r-^ m 

t=l t=l * t=l &m 2 T ~ t + 1 & 2 T 

for the second degree models, ignoring corrections of order a. An analogous result holds 
for the third degree model. Making a small angle approximation we arrive at 



n m 



t=i 



N tv* (38) 



where tv* is given in Tab. f|. Thus the product, being of order N 2t , balances the first 
factor which is of order a and hence N~ 2t . An orbit will be linearly stable if the remaining 
coefficient maximised over k (and denoted A*) is less than one. 

We can set /'(<&o) = by letting $ = for the second degree models and $ — — 1/2 
for the 3B model. For the 2A~ model we also have g'(0) = so the orbit remains 
superstable in the presence of perturbations at all wavenumbers k, and A* = 0. The B 
models have g'(0) = ±1, so that we have stability if A* = o*tv* < 1. This quantity is also 
given in Tab. |], and it shows that all orbits are linearly stable when /'($o) = 0. 

Note that the most stable state is when /'($) = belongs to a periodic orbit of the 
synchronised map (Q); this is not generally the same as the superstable point of this map, 
which is the solution of the equation (1 — a)/'($) + ag'($) = 0. 

Finally we make a connection to an observation made in Ref. 0. In section 6.6 of that 
work a scaling was observed numerically in the limit a — > 0, in particular the dynamical 
average 

v($) = (^/(<l> / )rf<^» , ) (39) 

Jo 

exhibited the following behaviour 

V{a) - 1/(0) = h{a)y/a (40) 

where the function h(a) satisfies (in the limit) 

h(N 2 a) = h{a) (41) 

that is, h(a) is log-periodic. Here we have found a family of stable orbits in which a 
differs from one to the next by a factor of iV 2 . The iV 2 comes from the derivative of the 
synchronised map at the points $ = ±1. The scaling of iV 2 would occur in other families 
of orbits (stable or unstable) containing $ = ±1, and hence possibly to all small values 
of a. 



3 Nonlinear stability 

Linear stability, considered in the previous section, implies that sufficiently small pertur- 
bations of a periodic state will approach that state in the future. However, this does not 
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imply that a generic initial state will approach the stable state. For a sufficiently large 
system, some of the initial conditions will be near points of the periodic state, however 
there is no guarantee that neighbouring maps will not strongly perturb these maps away 
from the periodic state. 

This section contains numerical results that use an initial state in which the $^0 
initial conditions are independently distributed uniformly from the interval [—1,1]. The 
spatial size of the system is 10 3 maps with periodic boundary conditions. After a time 
of 10 4 units the values of $ for each of the maps are plotted for many values of the 
parameter a. The results do not depend noticeably on the spatial extent, but do depend 
on the relaxation time (see below). See Fig. |4j. 

The boundaries of the linearly stable orbits of period 1 and 2 given in Tabs. § and |3] 
are shown as vertical lines; most of these boundaries are clearly bifurcations. The period 
1 orbits are both stable; the period 2 orbits while having a clearly observable effect on 
the dynamics, do not appear to be the entire attractor; the high period orbits discussed 
in Sec. ^]3| appear completely absent, although they are visible if the initial conditions 
are chosen closer to the periodic state. 

The most interesting question is what has happened to the period 2 orbits. Looking at 
the 2B~ in figure |], the period 2 orbit bounds a discrete set of values for 0.3675 < a < 0.48 
and a continuous set of values for 0.48 < a < 0.5556. In both cases a spatial slice of 
the solution reveals intermittent switching between the two values, as shown in Fig. [|. 
The discrete case contains dynamically "frozen" domains which remain for all time. The 
continuous case relaxes by increasing the size of the domains, but does so exponentially 
slowly. 

In the discrete case, there is no relaxation, so the final state (and consequently any 
property of it, such as the average V in Eq, ( ]3"5| ) depends strongly on the choice of ini- 
tial conditions; a different distribution function for the initial conditions would lead to 
a different distribution of domains. Averages can be estimated given a knowledge of 
the relative proportions of upper domains, lower domains and boundary regions. Other 
properties such as spatial correlation functions require more information about the dis- 
tribution of domains. 

The continuous case is numerically ambiguous due to the exponentially long time 
scales. It may may relax eventually to the stable periodic orbit, so its infinite time 
properties can be calculated directly from this orbit. 

As the relaxation time increases, the transition between the discrete and continuous 
regimes (ie a = 0.48 at a time of 10 4 ) is observed to move towards larger a, so it is in fact 
likely that at sufficiently long times, the final destination is a state with frozen disorder 
for all a corresponding to this stable periodic orbit. For example, the lower right plot in 
Fig. [5] shows that at a = 0.5 and a time of 10 8 the dynamics seems to have frozen. 

However, even if as suggested in Ref. the coupled map lattice describes something 
going on at the Planck scale (10~ 35 seconds), the transient relaxation may relate to 
measurable time scales. Note that time in the coupled map lattice is not physical time 
but Parisi-Wu fictitious time which supposedly is taken to infinity. The question is 
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a 

Figure 4: Bifurcation diagrams for each of the models. The left column from top are 
2A, 2A~ and 3A respectively; the right column is the same with A replaced by B. The 
vertical lines indicate boundaries of linearly stable orbits from Tabs. Q [| 
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Figure 5: Spacelike slices of the period 2 orbit of the 25„ map, showing frozen disorder 
(a = 0.4, left) and exponentially slow relaxation (a = 0.5, right). The relaxation times 
are 10 4 (top), 10 6 (middle) and 10 s (bottom) 
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whether this indeed happens or rather it remains exponentially large but finite. 

The other backward coupled (B) maps and also 2A_ have the same appearance of 
multiple "feathered" attracting points, which apparently give way to a continuous distri- 
bution as a is varied, specifically for IB and a > 0.6, for 2A~ and a pa 0.3 or a pa 0.7, and 
for 3B and a > 0.8. It seems likely that this phenomenon of exponentially slow dynamics 
and strong dependence on initial conditions holds for all of these regimes also. 

The other general feature, noted in Ref. f|, is that for small a, the advanced (A) and 
backward (B) coupled maps lead to very similar behaviour. 

Comparing Fig. |] here with the results in Ref. |J, the following points (from which 
the standard model parameters are calculated) are in regions with apparently stable 
attractors: 

First, the interaction energy zeros: ajj 3 ^ ~ 0.35 lies in the domain of stability of 
a period two orbit, however the numerical results suggest chaos here (for our uniformly 
distributed initial conditions), albeit strongly focussed on a small range of $. af A) pa 0.12 
lies just below the lower boundary of the stable window (of spatial period 2; not discussed 
above). ^ ~ 0.18 is close to the boundary of stability of the period 2 orbit at 1/6. 
0-2 ? ' A = 1/2 and a^ B ' = 1 both lie in the stable (period one) regions as noted in Ref. ||. 

Second, the self energy extrema: a' 3 / 4 pa 0.17 and a' 3 s A pa 0.23 both lie in a complicated 
stable region, as do a' 3B pa 0.19 and a'g 3B ^~ ' 29 , the latter of which is corresponds to a 
stable period two orbit. a'jj 2 ^ pa 0.18 lies close to the upper boundary of the stable 
window mentioned in the previous paragraph. a'^ B ^ ~ 0.22 lies in a stable period 2 
orbit, although numerical results suggest chaos in its vicinity. 

Thus, these values (and hence standard model parameters according to Ref. ||) relate 
to many different types of dynamics, from fully chaotic (ie not in the lists above) to a 
stable fixed point. Only the latter is easily obtained analytically. 



4 Conclusion 

Coupled Tchebyscheff maps provide excellent scope for both analytic and numerical stud- 
ies of stable synchronised states. Periodic states of (temporal) length one and two can 
be found analytically, and are found to have linear stability for a number of models in 
various ranges of the parameter a. Numerical observations confirm transitions at values 
predicted by the theory. 

There is also a family of orbits of increasing length that can be handled analytically in 
the long period limit, which demonstrates the existence of local stability, even arbitrarily 
close to the uncoupled fully chaotic limit a = 0. 

Numerically, the period 2 orbit in the 2B~ model is an attractor where the two values 
of $ occur in different domains at one time slice, exchanging their values at the next time 
slice. The size and distribution of these domains depends sensitively on the distribution of 
initial conditions. Furthermore, some values of a require exponentially long times before 
the final state is reached. Similar structures appear in the 2A~ , 2B and 3B models. 
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How do these results affect Ref. |J? In general we might expect the results to depend 
on the choice of initial conditions if a stable periodic state exists (true for many of the 
values) or indeed any coexistence of attractors (beyond the scope of this paper) , and also 
be affected by exponentially slow dynamics in one of the "feather" regions (apparently 
not relevant at the values given in Ref. 0). 

In addition, the fact that many if not most of the values of a correspond to stability, 
or at least not full chaos, undermines the assertion that it is the strong chaotic properties 
of the Tchebyscheff maps that is responsible (via stochastic quantisation) for quantum 
mechanics. The coupled Tchebyscheff maps are often far from chaotic, and some expla- 
nation must be given as to why many stable parameter values are important in particle 
physics yet the quantisation mechanism requires strong instability. 

The results presented here are only a brief sketch of the diffusively coupled Tchebyscheff 
map lattices, however it has been sufficient to find analytic results for classes of periodic 
states and demonstrate some of the richness of spatially extended dynamics including 
bifurcations evident in Figs. [I] and |], and exponentially slow dynamics in Fig. |5]. The 
combination of analytic tractability and a great variety of dynamical behaviour makes 
the coupled Tchebyscheff maps good candidates for prototypes of spatiotemporal chaos, 
much as the (closely related) logistic map is a prototype of low dimensional chaos. 
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